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The structural correlation functions of a weakly disordered Abrikosov lattice are calculated in 
a functional RG-expansion in d = 4 — e dimensions. It is shown, that in the asymptotic limit 
the Abrikosov lattice exhibits still quasi-long-range translational order described by a nonuniversal 
exponent rjc which depends on the ratio of the renormalized elastic constants k = cm/cu of the flux 
line (FL) lattice. Our calculations clearly demonstrate three distinct scaling regimes corresponding 
to the Larkin, the random manifold and the asymptotic Bragg-glass regime. On a wide range of 
intermediate length scales the FL displacement correlation function increases as a power law with 
twice the manifold roughness exponent Crm(k), which is also nonuniversal. Correlation functions 
in the asymptotic regime are calculated in their full anisotropic dependencies and various order 
parameters are examined. Our results, in particular the K-dependency of the exponents, are in 
variance with those of the variational treatment with replica symmetry breaking which allows in 
principle an experimental discrimination between the two approaches. 



PACS numbers: 74.60. Ge, 05.20.-y 
I. INTRODUCTION 

The discovery of -high- Tc superconductors in 1986 by 
Bednorz and Miilleiul has led to a strongly renewed inter- 
est in the theory of superconductivity, both in the expla- 
nation of its microscopic originH and in the application 
of the phenomenological theory on the determination of 
the phase diagram, of dissipation effects due to flux creep 
and related phenomenaoQ. 

Thermal fluctuations turned out to lead to a phase dia- 
gram drastically different from the mean-field prediction 
even in pure systems. This is due to the elevated transi- 
tion temperatures as well as due to the pronounced layer 
structure of the high- Tc materials. A liquid phase is now 
expected to separate the flux repulsing Meissner from 
the Abrikosov phase. In the latter the magnetic induc- 
tion B enters the material in the form of quantized flux 
lines (FLs) which still form a triangular lattice. How- 
ever, the upper critical field (where superconductiv- 
ity disappears) is now strongly reduced with respect to 
its mean-field value and can be understood as resulting 
from a melting of the flux line latticeB. In fact there 
are several experiments on high-Tc materials providing 
firm evidence that the Abrikosov lattice, is melted over a 
significant region of the phase diagramflfl. The field de- 
pendent melting line can be obtained from Lindeman's 
criterion in conjunction with a detailed knowledge of the 
elastic properties of the lattice. In general, these prop- 
erties have to be obtained bu. methods going beyond a 
simple continuum descriptiontJ. Moreover, at present, it 
is not clear, whether the transition to the normal phase 
at high fields happens in these materials via one or two 
transitions. 

It is well known that in addition to thermal fluctua- 
tions in type-II superconductors also the effect of disor- 
der has to be taken into account since FLs have to be 



pinned in order to prevent dissipation from their motion 
under the influence of an external current. Therefore, 
understanding the interaction of vortices with quenched 
randomness in form of these defects is of especial impor- 
tance. 

A natural question arising in this context concerns the 
influence of randomness on the translational order of tha 
lattice. This question was first considered by LarkinQ 
who found from perturbation theory that randomly dis- 
tributed pinning centers lead indeed to a destruction of 
the Abrikosov lattice. In particular, he obtained an ex- 
ponential decay of the correlations of the order param- 
eter for translational long range order (r) — e'*^"'^''^ 
on length scales larger than a disorder dependent Larkin 
length L^. Here G and u denote a reciprocal lattice vec- 
tor and the displacement field of the FL lattice, respec- 
tively, and r = (x, z) is a 3-dimensional position vector. 

This conclusion was in agreement with the more gen- 
eral observation of Imry and Ma, that quenched ran- 
domness coupled to a continuous symmetry order pa- 
rameter destroy true long range order in less than 
four dimensionaJ. In this context the Imry-Ma argu- 
ment played a similar role as the Mermin- Wagner the- 
orem for pure system in two dimensions. In the lat- 
ter case we know however, that the destruction of true 
long range order by thermal fluctuation contains still the 
possibility of topological order with an algebraic decay 
of correlationst3. The intriguing possibility of a weak 
disorder phase with a topological order, which distin- 
guishes it from the fully disordered phase at larger dis- 
order strength, remained as an open question, and was 
discussed intensively in recent works. ■— ■ 

As was first shown by Nattermanntil, in treating the 
interaction between FL lattice and disorder, it is cru- 
cial to keep the periodicity of this interaction under the 
transformation u ^ u + X, where X is a lattice vec- 



1 



tor of the Abrikosov lattice. This, symmetry, which is 
abandoned in peeturbation theorya and in the so-called 
manifold modelgl3, leads to a much slower, logarithmic 
increase of the elastic distortions, measured bjc _tiie_dis- 
placement correlation B{r) = ([u(r) - u(0)]2)EJ'yTa. If 
large dislocation loops can be neglected - as assumed in 
the above discussion - then the resulting phase is charac- 
terized by a structure factor with power-law singularities 
corresponding to a-ijuasi-long-range ordered flux phase, 
the "Bragg-glass"Efllj, with algebraic decaying pair cor- 
relation function Cdr) = (*g(i")*-g(0)). 

In three |-dipiensions there |— , is indeed 
strong theoretica£3~Ell and experimentalE3 support of a 
phase transition from a fully disordered phase, where the 
occurrence of unbounded dislocation loops leads to an in- 
stability of the Bragg-glass, to a dislocation-free genuine 
glass phase at weak disorder. 

The resulting power law decay of Cg (r) in this Bragg- 
glass is reminiscent of the situation in pure 2D-crystals 
where in the solid phase Cg(x_l) ^ |xj_|^'"=. This solid 
phase corresponds in fact to a line of critical points with 
?7g = TG'^{1 + K~^)/(47rcii) and k, = Cee/cn. The Cu 
represent the renormalized elastic constants which have 
a finite temperature dependent value. At the melting 
temperature Tm the exponent t^g reaches a non-universal 
value »«j = (6*0/471)^(1 — K^), which still depends on 

n{Tr,^)m. 

Based on this rough analogy one could expect a sim- 
ilar non-universal behaviour for the Bragg-glass phase, 
although it is dominated by randomness and, therefore, 
by a zero temperature fixed point. But, in addition to 
the relevance of metastable states, the accuracy of earlier 
renormalization group approaches and variational tech- 
niques was particularly hampered by the complex elastic 
properties of the Abrikpsow lattice. For instance, Gia- 
marchi and Le Doussallil'ta calculated Cg(x) using (i) 
a variational treatment for the triangular FL lattice and 
(ii) a functional renormalization group (FRG) ind = 4— e 
dimensions for a simplified model using a scalar displace- 
ment field u with isotropic elasticity only. In both cases 
they found Cg(x^,0) ~ |xx|^''"^ with a universal eyi^o- 
nent 7?Go = ^(4 - d) with A = I and A = « 1.1 

for the treatment (i) and (ii), respectively. Here Gq de- 
notes one of the smallest reciprocal lattice vectors with 
Goa = ^■nj^fZ, and a = (20o/'\/35)^^^ is the lattice spac- 
ing. To date, the question to which extent these results 
depend on the applied techniques and the simplifications 
of the considered models is not fully resolved. 

In this paper, we systematically study the structural 
properties of the triangular Abrikosov lattice at weak 
disorder using a FRG method. Contrary to previous 
approaches, we explicitely take into account the trian- 
gular symmetry of the lattice and all of its elastic modes. 
We derive functional recursion relations for the correla- 
tion function of the random potential, which, in combi- 
nation with a Fourier decomposition technique, allow us 
to extract detailed information about the collective wan- 
dering behaviour of the FLs. This transversal wandering 



can be characterized by the roughness exponent C, con- 
trolling the displacement correlations via B{y) ~ |rp^. 
Following the KG flow, three different scaling regimes 
can be clearly identified: On length scales smaller than 
the Larkin length, the FL are displaced by an amount 
smaller than the characteristic scale of the short range 
correlated random potential. Thus, Larkin's perturba- 
tion theory holds with Q — (4 — (i)/2. Beyond this regime, 
but for line displacements still smaller than their distance 
a > ^, one enters an intermediate regime with non-trivial 
roughness influenced by metastable states. Although this 
regime is similar to a single FL system, the complex elas- 
tic interactions of the lattice lead to a non-universal C 
depending on the ratio cge/cn as reported here for the 
first time. Finally, on asymptotic scales, the FL displace- 
ment becomes larger than a, leading to the logarithmic 
roughness responsible for the quasi-order of the Bragg- 
glass phase. Using our Fourier technique to analyse the 
FRG flow equations, we were able to determine from the 
fixed point value A* of the variance of the random po- 
tential the exponent t^g = A*(Ga)^ with sufficient accu- 
racy to rule out universality. Instead, rjc, depends also 
on the ratio cge/cn. The situation is therefore indeed 
qualitatively similar to that of 2D crystals at the melt- 
ing temperature as speculated above. With the ratio 
cee/cii depending in general on B and T, the observa- 
tion of a field-dependent 77G would yield the opportunity 
to judge the ptlidity of different approximation schemes 
under debatcEll. A brief description of our combination 
of FRG techniques and Jipmerical calculations and these 
results appeared earliernJ. 

The paper is organized as follows. First, we will intro- 
duce the model and examine its symmetries in Section 
2. The derivation of the FRG flow equations and its nu- 
merical treatment is outlined in Section 3. A detailed 
discussion of the different scaling regimes and its corre- 
lation functions is given in Section 4. In this Section we 
also study correlation functions of various order parame- 
ters to obtain information beyond the translational order 
parameter as a signature of residual order at weak dis- 
order. Experimental implications are presented briefly 
at the end of Section 4. Finally, in the last Section we 
summarize and discuss our central results. Lengthy com- 
putation is documented in the appendix. 

II. THE MODEL 

A. The Hamiltonian 

In order to examine the influence of disorder on the flux 
line lattice we employ the elastic description reviewed in 
Ref. |. The de grees of freedom are the positions of the 
flux lines {Rx(2)} in the [x, j/)-plane as a function of z, 
the direction of the external magnetic fleld, and X, the 
lattice vectors of the perfect triangular Abrikosov lattice. 
To measure the disorder induced distortion of the flux 
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lines, we introduce the displacements ux(^) = Rx('Z)— X 
from the perfect lattice. Such a labelling with the perfect 
lattice positions is possible provided no dislocations are 
present in the system. The stability against dislocations 
has to be assured a postcspiji-aihicl^^ a finite limit 
for the disorder strengtlTLMOaElM. 

In an elastic continuum description the displacements 
are extended smoothly to a continuous function ux(z) — > 
u(r) with r — (x, z). Syjametry reasons require three 
distinct elastic coefficientsEEl and the energy cost of a dis- 
tortion of the triangular line lattice is 

We; = ^ y" d^a^dz |cn (V^u)' + cge (V^ x u)' + 

+C44(V,uf}. (2.1) 

The elastic coefficients correspond to compression (cn), 
shear (cge) and tilt (044) of the lattice. The former two 
describe flux line interaction, while the tilt modulus con- 
tains both an interaction contribution and individual line 
tension. Such an elastic continuum description is valid 
as long as both displacements and their gradients vary 
slowly over lattice steps, which can be subsumed in the 
assumption of La ^ a. Here, La is the positional corre- 
lations length, defined as the distance of two FLs whose 
mean displacements vary by the order of one lattice spac- 
ing a. In our context this is a condition on the weakness 
of disorder. Below we show that it can be easily satisfied 
by realistic impurity concentrations. 

In wavenumber space, the elastic energy reads 

= i / (Sl'^^ + St'Pt) u(-q). 

BZ 

(2.2) 

The integration runs over the Brillouin zone (BZ) and 

are projectors onto the lon- 
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gitudinal and transversal modes, respectively, with prop- 
agators 



r'-l 2 



(2.3) 



This Fo urie r space formulation is more general than the 
model ( ^.l[ ) since it both allows for wavenumber depen- 
dent elastic moduli and docs not include the continuum 
limit as only wavevectors within the first Brillouin zone 
are included. The derivaiifm, of the elastic moduli from 
Ginzburg-Landau theoryEJ'cS indeed gives g-dependent 
elastic constants on scales smaller than the London pen- 
etration length A. However, this only leads to weakly 
renormalized constants on scales larger than A. For the 
case of weak disorder considered here, we focus on the 
behaviour on scales far beyond A and may thus content 
ourselves with the local version {2A). The coefficients 
Cii are understood to be the renormalized ones. In Eq. 



.2D we have formulated the model in general dimen- 
sions and extended the z-direction to a d—2 dimensional 
space. This will allow for an expansion around d = A in 
the renormalization procedure below. 

The interaction of flux lines with impurities competes 
with elasticity and tends to roughen the lattice. It is 
modelled by the coupling 



X 



= E /dz'*"'l^*.(X + u(X,z),z) 



of the flux line density 

PFLL(r) -^<5(x-X + u(X,z)) 



(2.4) 



(2.5) 



X 



to the disorder potential Vdis(r) with short range corre- 
lations on the scale ^ = max(^sc, £,dis), where £_sc, idis are 
the superconductor correlation length and the correlation 
length of disorder density fluctuations, respectively. The 
disorder potential is taken to be a random variable with 
Gaussian distribution. It is characterized by 



Vd^s{v) = , Vd^s{r) Vd^s[r') = A^(x - x') 5{z - z'), 

(2.6) 

where A ^ (x) is a delta- function smeared out over a region 
of size ^. 



The continuum limit of Eq. (2.4) has to be taken with 
some care. Slow variations of the displacements do not 
suffice to allow a straightforward continuum limit. The 
potential Vdis(r) varies rapidly over lattice steps since 
^ <C a. Therefore, in the Hamiltonian rewritten using 
Poisson's summation formula 



5]/(X)= /dV(x)[l 
X •' 



the terms corresponding to reciprocal lattice vectors G 7^ 
are important for tjicjiinning energy and have to be 
considered adequatelylHlo. 



dis 



Po 



Po 



Vdis (x -t- u (x, z) , z;) [1 + 



Vdjs(x, z)[l 



giG[x-u(x,^)]]^ ^2.7) 

G^O 



where po = B/$o is the average FL density. In the last 
equation, we have substituted 



X = X -t- u(x, z) 
u(x, z) = u(x + u(x, z), z) = u(x, z). 



(2.8) 



With u(x), the displacement field is now written as a 
function of the actual position x of the flux lines rather 
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than as a function of the perfect reference lattice coor- 
dinates X. Such a relabeUing is possible provided that 
|Vmq| < 1, which is assured by the above assumption of 
\WuJ\ i2 a/La < 1. 

The change of parameterization in Eq. (2.8) from 



{x, u} to {x, u} shall also be performed in the elastic 
energy via 



9x 



di 



(2.9) 



This gives one formally invariant term plus terms of 
higher order in derivatives of u. The latter are - due 
to the extra gradients - irrelevant in the RG to follow 
and therefore neglected in the first place. The same rea- 
soning holds for the extra gradient terms resulting from 
the change of measure in the partition function of the 
problem, where originally one sums over the configura- 
tions u(x). With this in mind, we drop from now on 
the tilde, keep the elastic lattice energy unchanged and 
derive a com pact form for the pinning energy. Retaining 
in Eq. (2/7) only the lowest order contributions in daUp 
and subtracting u-independent terms, we obtain 

Hms ^PO [ {- Vd^s{r)\7^u{v) + Vd^s{r) ^ e^G[x-u(r)]| 

(2.10) 



The first term couples the divergence of the displacement 
field to the disorder potential. It can be shown to lose 
against elastic energy in the effort to roughen the lattice 
above two dimensions by a simple scaling argument. One 
assumes the system to be rough and the displacement to 
vary with the scale as u ~ L^, C < 1- The elastic energy 
then scales as L''~2+2C g^j^j th e pin ning energy from the 
first and second term in Eq. (2.10) scale as ^(^^-2+20/2 
and L''/^, respectively. Therefore, for d > 2 the first term 
can be neglected against the elastic energy, whereas all 
terms of the second term are relevant for d < 4. The 
relevanc e of t hese infinitely many exponential operators 
in Eq. ( 2.10 ) makes a functional renormalization group 
inevitable. The preceding discussion shows that the 'cor- 
rections' to the naive continuum limit are dominant and 
retaining the lattice structure is crucial (see also Section 

iil(c)). 

Conse quen tly, we confine our analysis to the last term 
in Eq. (2.10). The slowly varying displacement field in 
the harmonics rather than in the disorder potential itself 
better shows which terms are relevant. This is made use 
of when we average over the disorder potential, which 
is done by the standard replica method. The averaging 
process introduces an interaction between different repli- 
cas, which is calculated in detail in Appendix ^ The 
resulting replica Hamiltonian, which will be the starting 
point for the renormalization group (RG) below, reads 



1 " r 
1 " r 

E d'rR{u^{r)-u\r)), (2.11) 



with the disorder correlation function defined by 

3-*^". (2.12) 



i?(u)^pg^ A(G) 



The random potential correlator A^(x) [with Fourier 
transform A(k), see Eq. (2.6)] is given in terms of phys- 
ical quantities by 



(2.13) 



where fpi„ is the mean individual impurity pinning force, 
TT-imp the impurity density and g(x/^) a function of am- 
plitude 1 and range ^, cf. Ref. This allows to model 
the bare, unrcnormalized correlator as 

i?o(u) = /2^„n.„pe'|^ E - GOe'°". (2.14) 



G 



The second derivatives Rap = du^du^R of i?(u) at the 
origin will be of central interest below. Using Eq. ( 2.14 ) 
its unrcnormalized bare value can be approximated by 
i?,40) = i?,,(0) ^ 102/p'.„^™p^'S/*o (1 + 0{S,/a)). 



B. Symmetries 

Next we want to examine the symmetries of the system 
in Eq. (2.11) and their constraints on the structural order 
of the flux line lattice. 

(a) TLei is isotropic in the (x, j/)-plane. For the special 
choice of cii = cge one can make TLei isotropic even 
in (x, z) by repealing z z^Jc^^Jcxi. Previous 
RG approaches^ were limited to this very special 
choice ignoring anisotropy. 

(b) The elastic energy Tiei is invariant under si- 
multaneous rotation of displacements by 90 de- 
grees (u^(Y),Uy(Y)) {uy{r),-Ux{r)) and k = 
cee/cii K^^. Divergenceless shear strains are ro- 
tated into curlfree compressional strains, compen- 
sated for by the exchange of the respective moduli. 
In wavenumber space the longitudinal modes are 
exchanged with the transversal ones. 

(c) The disorder correlation function R{u) (see Eq. 
( ^.12 )) has the full symmetry of the triangular lat- 
tice, i.e., it shows invariance under translations of 
u by lattice vectors and rotations of u by integer 
multiples of 60 degrees. The former is obvious from 
i?(u) being a Fourier sum while the point group 
symmetry of i?(u) can be seen easily if one applies 
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a rotation D, that maps the lattice onto itself, and 
uses the isotropy of disorder 



i?(Du)=p2^A(G)^ 



-iGDu 



G 

DG' 

= pIJ2 AiG') 

G' 



-iG'u 



-iG'u 



While disorder is assumed to be homogenously and 
isotropically distributed, the reference state of the 
FL positions shows the discrete triangular sym- 
metry. Together, for the displacements measured 
from the lattice, this gives a disorder correlator 
with triangular symmetry. Due to the discrete 
6-fold rotational invariance the function R{u) be- 
comes isotropic at the origin, i.e., Rxx{0) — Ryy{0), 

Rxy{0) = 0. 

(d) The effective replica pinning energy is invariant un- 
der 



u''(r) u"(r) +f(r) 



(2.15) 



for an arbitrary (yet constant in the replica index 
a) function f(r). This is a consequence of local- 
ity of our compact form for Tidis and thus only an 
approximate symmetry. As can be seen from the 
detailed analysis in appendix ^ locality grows if 
A^(x) gets sharper. Upon renormalization A^(x) 
becomes more and more delta-function like. This 
is not surprising as the short scales, whose coupling 
by a finite correlation gives the non-locality, are in- 
tegrated o ut. O n larger scales we may thus take 
symmetry ( 2.15| ) as given. 



These symmetries allow for important conclusions, which 
simplify the calculations to follow: 

(i) Symmetry (d) grants that the elastic constants in 
Tiei are renormalized only trivially by rescaling in 
a RG procedure. This is a well known property 
of systems consisting of an elastic term diagonal in 
replicas and an interaction term, that depends only 
locally on differences of repfeiajiclds. It is often ref- 
ered to as 'tilt symmetry'c3E3 although here it is 
generalized to tilt, shear and compression moduli. 
As the symmetry is fulfilled in the present case only 
for larger scales, elastic constants will be renormal- 
ized weakly on small scales. We start our descrip- 
tion there and the effective constants shall for con- 
venience again be denoted by cn, C44 and cqq. 

(ii) The prime measure for residual translational order 
in the system will be the mean squared relative 
displacements 



Their scaling Bapir) ~ jrp'' defines the roughness 
exponent C of the lattice. From (a) and (c) follow 
the isotropic relation Bxx{x,y, z) = Byy{y, —x, z) 
and analogues. 

(iii) Symmetries (b) and (c) give B^^ (r) = B^^ ^ (r) and 

^xyi^) — ~^xy (j")! which are relations between 
correlation functions of different flux line lattices 
which are related by an inverted ratio k = cge/cn 
of elastic constants as indicated by the superscript 
on B'^p{v). 



III. FUNCTIONAL RENORMALIZATION 
GROUP 

As was motivated in the last Section, in d < 4 di- 
mensions we have to deal with infinitely many relevant 
operators. Therefore, we employ a functional renormal- 
ization group (FRG) method to treat the interaction po- 
tential i?(u) between different replicas. In the past, this 
technique had been successfully arai 
tic objects in random environmentstiJ' 
method was usedlp describe lattices of elastic objects like 
that of flux linesEj. However, the model studied by FRG 
for the latter case has two important drawbacks: The ef- 
fect of the triangular lattice is neglected in both (i) the 
elasticity, that is anisotropic for all physical FL lattices 
and (ii) the disorder correlator. The latter reflects the 
lattice symmetry and any RG flow will have to preserve 
it. In the following we develop a FRG approach to take 
into account both effects, which result in new physical 
behaviour. 



;d to single elas- 
I. Later the same 



A. RG equations 

We use a standard hard-cutoff RG by integrating out 
the displacement field u(q) with wavevectors q in an in- 
finitesimal momentum shell below the cutoff A > |q| > 
A/6 = Ae"*^'. The contributions to the RG equations 
from rescaling according to 



Y — y' h 
u(r) = u'{Y')b^ 



are given by 



f L=(2-rf-2C)T 
dR 



(3.1) 



(3.2) 



B^p{y) = (K(r) - u„(0)][u^(r) - ^^(0)]). 



^ Lc = (4 - - 4C)i?(u) + C9ai?(u)uo. (3.3) 

Since the order of the flux line lattice is expected to be 
dominated by the random potential, we have chosen to 
rescale the temperature in order to organize the RG anal- 
ysis of the expected T — Q fixed point. The corresponding 
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Eq. (3.2) is exact, i.e., there will be no feedback from the 
random potential to the elastic constants as mentionend 
above. From the second term on the rhs of Eq. ( |3.3D it 
is obvious that the periodicity of i?(u) is not consistent 
with a finite roughness exponent C. Therefore, we have 
to choose ^ = to obtain a periodic fixed point function 
i?*(u). This choice of C automatically implies smaller 
than power-law roughness on asymptotic scales if a fixed 
point is found. Still, choosing C to be zero is just a mat- 
ter of the fixed point analysis and does not forbid finite 
values for ( on smaller length scales. 

To obtain a fixed point for i?(u) to order e = 4 — d, we 
have to calculate a ll te rms of second order in i?(u) which 
contribute to Eq. (3.3). To do so, we trace over the fast 
modes of u(r). The feedback to the disorder energy can 
be written in a cumulant expansion in Tidis ■ Taking into 
account the anisotropy of the elastic kernel, we obtain 
the RG equation 



dR{u) 



diR{u) = ei?(u) + ( -d^d^Riu)dpdsRiu)- 



dad-yR{u)dpdsR{0) ) M"^^^'' (3.4) 



with 



= dl 



r gMn)G^si-q) (3.5) 



with Gap = GL^ap + GrPap- Hcrc denotes integra- 
tion over the shell A > |q| > Ae"'^' to order dl. The 
symmetric matrix M with the greek double index ordered 
according to {xx, xy,yx, yy) is computed in appendix ^ 
as 

3/1 + 2/2 /1+6/2 

1 I /i - 2/2 /i - 2/2 
/i - 2/2 h ~ 2/2 

/1 + 6/2 3/1 + 2/2 



(3.6) 



with integrals 



h = dr^ 



r + Gi) = 

•la 



h = dl- 



GtGl 



1 1 + K 

Stt^ C44C66 

1 In K 
87r2c44Cii K — 1 



(3.7) 



In the spirit of a consistent e-expansion, the integrals are 
evaluated in 4D. 

With the shorthand notations —A = Rxx{0) = 
Ryy{0), 6 = 1 — 2/2 //i, and with the rescaling 



R = R— — . 

h 



(3.8) 



which makes the parameter A dimensionless, the RG 
equation finally becomes 



dl 
-,2 



eR{u) 



+ y {7?L(u) + Rlyiu) + 2Rly{u) + 2A[Rx4n) 
+Ryy{u)]~-^-{[RUu) - Rlyiu)]' + iRlyin))}. (3.9) 



It is easily checked that the RG flow given by Eq. (|3.9| ) 
preserves the symmetries of the correlator R{u). The RG 
flow depends on the ratio k = cqq/ch via the anisotropy 
parameter 



5=1-2- 



In K 



K — K 



(3.10) 



which varies between and 1. This dependence on the 
elastic constants is inherently related to the anisotropic 
elasticity and cannot be eliminated by further rescaling 
as it is the case for the more isotropic idealizations of 
flux line lattices studied in Ref. |l^. Below it will be 
shown that the nonuniversality of the coefficients in the 
RG Eq. (3.9) carries through to nonuniversal exponents 
for the displacement correlations. To compare with for- 
mer studies, we consider two limiting cases of Eq. (3.9): 

(i) cii = C66 corresponding to (5 = 0. This is the isotropic 
case which former studies were restricted to, mainly for 
technical reasons, since this extreme limit is not realized 
in isotropic superconductors at low temperatures. 

(ii) cii ^ C66 corresponding to (5 = 1. This case is real- 
ized physically near the upper critical field Hc2 where the 
lattice becomes very soft with respect to shear stress. In 
the very extreme limit cn — > c» only the transverse prop- 
agator is effective in the Hamiltonian and Eq. ( |3.9D can 
be rescaled as to be free of the remaining elastic constant 

C66- 



B. Numerical solution 

The quantity which will determine the effective propa- 
gator of the displacement field is the renormalized disor- 
der strength A/ , considered as a function of the logarith- 
mic l engt h scale or RG parameter I. Since the RG flow of 
Eq. ( 3.9 ) cannot be reduced to that of A/ only, we have to 



solve the full partial differential equation for the function 
Ri{vl\ starting from its bare value i?o(u) given by Eq. 
(^.14 ) . Even for the fixed point condition dR{\i) /dl = 
an analytical solution is not obvious. Therefore, wc treat 
the equation of flow numerically. Rather than solving the 
fixed point equation, we obtain the fixed point by inte- 
grating Eq. ( [3.S| ) numerically from the starting function 
i?o(u). This is necessary both to exclude non-physical 
fixed point solutions, which are not connected by an RG 
flow to i?o(u) and to get additional information about 
the b ehaviour on intermediate length scales, see Section 
IVA. 



As mentioned above, the lattice symmetry of the initial 
function i?o(u) is preserved by the RG flow. Technically, 
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this allows easily for a numerically stable and effective 
iteration of the RG flow to large I by rewriting i?(u) as 
the Fourier series 



R{u) = ^ i?G cos(Gu) 



(3.11) 



G 



with identical coefhcients for reciprocal lattice vectors re- 
lated by a pointgroup symmetry. The partial differential 
equation for i?(u) becomes with this ansatz an infinite 
set of coupled ordinary differential equations for the co- 
efficients {Rg}. Of course, the numerical integration has 
been restricted to a finite set of coefficients by choos- 
ing a cutoff Gc so that |G| < Gc- The accuracy of the 
computation depends on Gc- To obtain numerical re- 
sults for Gc — > oo, we employ a routine typically used 
for finite size scaling of numerical data to get sufficient 
precision for dxxR*iO) — dyyR*{0) = — A*, the quantity 
entering further calculations. How reliable the numerical 
results indeed are, can be estimated from the exact rela- 
tion for the fixed point correlator eR*(p) + a'^{A*f = 
that stems from Eq. (3.9), see Figure H. The accuracy of 
4 significant digits of A* can be assured. 

The variable parameters of the initial function Rq{u) 
are the number of non-zero coefficients Rg , which is fixed 



by the value of ^ in Eq. (2.14), and the magnitude of the 
coefficients Rg- Notice that all Rg are driven to non- 
zero values by the RG flow even if their bare values were 
choosen to be zero at the beginning. Whereas the choice 
of ^ in relation to a is important for the existence of an 
intermediate length scale regime with a non-trivial finite 
C (see Section IV A), the magnitude of the Rg deter- 
mines just the crossover length scale to the asymptotic 
fixed point function i?*(u), which itself is universal. This 
universal behaviour has been confirmed numerically for a 
very wide range of starting values. One can write the con- 
dition for finding convergence to the fixed point as a lower 
limit for the positional correlation length La > lO^A^^. 
A is the small lengthscale cutoff and of order A, the Lon- 
don penetration length. This condition can now be com- 
pared to 

which is the weak disorder condition necessary for the 
system to be stable against the formation of dislocation 
loopsEj. C is a constant of order one and Oc the core 
radius of a dislocation. One recognizes that the bassin 
of attraction of the fixed point covers well the range of 
validity of our dislocation-free description. 

In Figures H, ^ the fixed point function and its sec- 
ond radial derivative are shown. Fig. ^ shows clearly 
the existence of a non-analytic cusp at the lattice sites. 
Dependency of i?*(u) on the anisotropy is obvious from 
the numerical results for A* as shown in Fig. |3[ 



IV. CORRELATION FUNCTIONS 



-R'(0)-a'A'' 
, -R"(0) 



12 56 192 360 

number of non-equivalent coefficients 

FIG. 1. Numerical accuracy as a function of non-equivalent 
Fourier coefficients, i.e. coefficients that are not related by 
a point group symmetry. Our extrapolation method yields 
for the test quantity (e^*(0) + a^{A*f)/eR*{0) = 6 • 10"'^ 
instead of the theoretically expected value of zero. This gives 
the precision for A* mentioned in the text. 




FIG. 2. Fixed point correlator R*{u)- The hexagon repre- 
sents a Wigner-Seitz cell of the triangular lattice. 



A'/e 




FIG. 3. Variation of A* = -_RJ^(0) = -Rt,y{0) with 
anisotropy 5 = 1— 21nK/(K — k~^). 
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A. Length crossover effects 

With the solution of the RG equation at hand, we cal- 
culate the disorder and thermal averaged squared rela- 
tive displacement as a measure of translational order. It 
is related to the basic correlator of the replica fields 

K{ci)uli<i')) = TG^J^^ia) (2vr)'^<5(q + q') 

by 

Ba/3(r-r') 



ln(L|A) 



ln(L„A) 



{[uair) - u„(r')][u/3(r) - u,}(r')]) 
21im fe^(q) (l-cosq(r- r')) . (4.1) 

n^O J 



If we denote the correlator for a system with a renor- 
malized Hamiltonian Tif by Qi and if we allow tempera- 
ture to flow according to Eq. (3^), the scaling relation 
Tg{q) = Tigi{bci)b'^^+'^ holds for aU q with |q| < A/5. To 
obtain to order e the exact propagator t/(q) for fixed q, 
we choose h — A/|q|. This allows to employ a harmonic 
approximation to calculate Qi since the coupling of modes 
with momenta between |q| and A has been taken into ac- 
count already by the renormalization of 7i". Thus we 
obtain the propagator 



rf-4 

-1 A, 




(ei(q)P'^ 



a(3) 



-T (gL(q)P 



a/3 



(4.2) 



with I = InA/g. The term T is the mere thermal 
propagator. It does not effect roughness in d > 2 below 
the melting transition and is thus negligible against the 
disorder term ~ A; which is more strongly divergent for 
|q|-0. 

Before we discuss the asymptotic displacement corre- 
lations in the next Section, we want to exploit the be- 
haviour of Bap{v) on intermediate length scales. Ac- 
cording to the definition of the roughness exponent C 
by B{v) ^ |rp^ the scaling of the renormalized disor- 
der strength A/ determines C^. It is easily observed from 
Eq. (4.2) that necessarily In A; ~ 2C,l before it reaches its 
fixed point A* corresponding to C = 0. The RG flow of 
In A; as a function of I is plotted in Fig. ^. Two qualita- 
tive different types of behaviour can be observed, depend- 
ing on the bare form of the correlator i?o(u): (i) If a = ^ 
only the lowest order harmonics of Ro{m) are non-zero 
and two different scaling regimes with a sharp crossover 
emerge. The first one is called random force (RF) regime 
since C, = Crf = e/2 as predicted by Larkin's perturba- 
tive approach for a random force model3, which becomes 
applicable if B{y) < The last scaling regime is the 
asymptotic (A) one with logarithmic roughness, (ii) If 
a > ^ also higher order Fourier coefhcients are non-zero 
and a new scaling regime with a non-trivial value of C 



RF 

-<E >- 



RM I asymptotic regime 

^ 3:-l<^ 



A, 
10-= 

10-^ 
10-= 



5 10 15 20 

I = In (LA) 

FIG. 4. Scaling behaviour of the renormalized disorder 
strength A; in the three different length scale regimes (solid 
curve). The random manifold (RM) regime disappears 
(dashed curve) if the bare iio(u) has non- vanishing Fourier 
coefficients for the lowest order hexagon only, corresponding 
to the case a = ^ [see Eq. (2.14)]. 





FIG. 5. The roughness exponent ("rm is nonuniversal in the 
RM regime. 



appears between the RF and A regime. It is called ran- 
dom manifold (RM) regime since here the roughness is 
expected to be determined by the wandering of lines that 
do not yet compete with each other for disorder energy 



minima, i.e. 



B{v) 



< 



But metastability is already 



important on these length scales leading to a new rough- 
ness exponent Crm- One of the central results of our 
analysis is that Crm is nonuniversal. It depends on the 
ratio K of elastic constants and varies between 0.1737e 
and 0.1763e as shown in Fig. ||. The only result for Crm, 
which has been available so far, is based on a Flory type 
argument for the flux line lattice in that regimeOt^. It is 
given by Cam — — e/6 — 0.167e {N is the number 
of components of the displacement field, here N = 2.) 
and supposed to be a lower bound, in agreement with 
our findings. We expect that the range of variation of 
the nonuniversal Crm is larger in real three dimensions 
than an epsilon expansion to first order can reveal. This 
is because in the propagator the effect of variing cn, cgg 
is suppressed in an expansion around 4D by the more 
heavily weighted C44q^-term. 

Next we calculate the crossover length scales. These 
are denoted by and La for the crossover RF - RM and 
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RM - A, respectively. can be obtained from Larkin's 
perturbative analysis, which yields for the displacement 
correlations 



B(r) = ([u(x,z)-u(0,0)]2) 



$0 



1/2 

C44 



X [cfj\^- + z^f/^ + c-J/\.' + z-f'^] . (4.3) 

Here we have introduced the rescaled z-coordinates z; — 
z-\/cii/c44, Zt = Zy/c^oJciA- The anisotropic crossover 
or Larkin length scale is determined by the conditions 
S(0, z = L|) ~ ^2 and B(|x| = L^, 0) ~ , thus giving 



LI 



7r$o 



C11C44C66 



BfpinlT-imp Cll + C66 

7r$o C4f (ciicee)^/^ 



3/2 3/2 



-11 



(4.4) 



66 



Based on the condition for the mean displacements of 
lines at the crossover scales, the second crossover length 
scales L^, should be related to the Larkin length 



scales by La 



(^,x) 



(a/^)i/CRMi 



To check this rela- 



tion we have determined the crossover scales numerically 
from the RG flow of A; for different ratios a/^. The es- 
timates for C.RM thus obtained are consistent with the 
exact values shown in Fig. |5[ Of course, due to the fi- 
nite extent of the crossover regions, errors are much too 
large as to take this as a measurement of C.rm] however, 
it clearly confirms the physical picture of the crossover 
lengthscales. Turning to the RG flow of the whole func- 
tion Ri (u) , one expects that the appearance of metasta- 
bility on the Larkin scale should be refiected in a change 
of the functional form of i?;(u), too. Indeed, the exis- 
tence of a cusp in the RM and A regime can be observed 
from the sequence of Figures ^, 0, |g and ^, which show 
the second radial derivative of Ri (u) in the different scal- 
ing regimes, see Fig. IJ. From this observation we expect 
that the fourth derivatives dt^R{0) = d^^R{0) = M^^O) 
diverge at the scale given by the Larkin length. The 
scale Ic — In(icA) at which t he c usp appears can be 
extracted from the RG Eq. ( |3.9[ ) as pointed out in 
Ref. p3 for a different RG equation. Using the relations 
d^RiO) = 5*^(0) = 39292^(0) amongst derivatives of 
the correlator at the lattice sites, we get as equation of 
flow for the fourth derivative (0) at the lattice sites 



9,i?(4)(0) = ei?(4)(0) + —{10-S) R(^\o) 



(4.5) 



This equation can be easily integrated and yields in 
the weak disorder limit for the length scale where first 
Lt^iO) = 00 the resuh 



L. 



1 



a2A(10-<5)^W(o)- 



(4.6) 



Using relation (2.14), (p.^) the bare value of rI^^ can be 



expressed in physical quantities and we get for f <C a 



1 



1 



C44C11C66 



2po/, 



Cll + C66 



(4.7) 



Com parison of this result to the Larkin length of Eq. 
(4.4) shows that the correlator i?;(u) becomes indeed 
non-analytic at the crossover to the RM regime. 




FIG. 6. Second radial derivative d^^^R{u) at the initial po- 
sition a (Z = 0) shown in Fig. W. 




FIG. 7. Random force regime: d^^^R{u) at position bin 
Fig. ^ obtained by integrating numerically the RG Eq. (3.9). 



B. Asymptotic displacement correlations 

We now study the correlations in regime A in detail, 
starting with the displacement correlations. Beyond the 
second crossover scale La the larger-g modes governing 
both the RF and RM regime become unimportant and 
we may benefit from the fact that A; — > A* in the asymp- 
totic regime. Therefore, the only relevant disorder part of 
the correlator in Eq. (|!|) behaves like T^°^(q) ~ A*^-'^ 
leading to logarithmic roughness. It is interesting to 
note that the amplitude of the displacement correlations 
Bapir) does not depend on the disorder strength, but is 
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and thus 



FIG. 8. Random manifold regime: 9j^^|-R(u) at position c 
in Fig. a 




FIG. 9. Asymptotic fixed point; 9^^|i?(u) at position d in 
Fig. |. 



proportional to the fixed point value A* (k) , whose depen- 
dency on K is shown in Fig. ^. The correlation function 
Bapi'^) can now be calculated with expl icit co nsid eration 
of the anisotropic elasticity using Eqs. (4.1), ([4.2|). With 



the abbreviations zi = cn/ caA: zt = zy^cee/cAA and 
h{t) = t~'^ ln(l + t^) we find in regime A 



Byy{r) = 



A*(K)a^ 
1 + K 

9 ? 



1 + K 

9 2 

X ~ y 



In 



LI 



K In 



1 — K — h 



Zt 



K h 



In 



K In 



1 — K — h 



Zt 



K h 



2A*(K)a2 xy 

1 + K 



K — 1 — K 



Zt 



(4.8) 



([u(r)-u(0)]2) =i?,,(r 
2A*(K)a2 



1 



In 



£2 



K In 



+ zf 
LI 



(4.9) 



The details of the Fourier transformation are given in ap- 
pendix |C| There the calculation retaining all unbounded 
terms for large |x|, z is performed in 4D since A* is ©(e) 
For the reason mentioned in Section 
be obtained from B^xi^) by k — > 



II B 



(b) Byy{v) can 
The relative 



line displacement Bapir) grows only logarithmically due 
to the restriction of available configurations by adjacent 
lines. The 'geometric' coefficients of the logarithms re- 
fiect the elastic anisotropy. The result reminds of the 
findings for the analogous 2D flux line lattice by Carpen- 
tier and Le Doussal, where there is also a nonuniversal 
prefactor of logarithmic roughness that depends on the 
elastic details of the systemEj. 

To compare our results witli-former results by Gia- 
marchi and Le Doussal (GD)Eil^E3 in 3D we consider two 
limiting cases. First, in the limit z ^ the relative dis- 
placement reduces to 



— K — 1 



S,, = A*(K)a^|ln^^— 
2 K - 1 



Bxy = 2A*(K)a' 



X^ K + 1 



This result is identical in form with the result of GD ob- 
tained from a variational ansatz treatment for the trian- 
gular lattice. However, this ansatz yields a value for A* 
which is independent of n and differs from ours by about 
15%. The origin of this deviation is the inexact treatment 
of fluctuations by the variational method. The renormal- 
ization group approach of GD is restricted to an idealized 
scalar field model of a flux line lattice with isotropic elas- 
ticity. To determine the influence of the triangular lattice 
symmetry on A* compared to the scalar model, we con- 
sider now the isotropic case Cn = Cgg corresponding to 
K = 1. With the rescaled coordinate r' = (x, ch/caaz) 
we obtain 



Bxx ^ Byy ^ 2K* {1)0" \n{\r'\/ La), 



Bxy = 



with A*(l) = 0.0217e. This has to be compared to the 
result A*(l) = e/36 ~ 0.028e for the scalar model. 

C. Order parameters 



Now we study the translational order parameter 
^'g(i") = e'*^"*-""-* to measure the remaining translational 
order in the impure system. The pair correlation function 
is 



Cg{y) = (*G(r)*£;(0)) = (e<G(u(r)-u(o))), (4.10) 
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It is a preferred measure, as its scaling behaviour deter- 
mines the intensity of the reflection pattern obtained in 
neutron scattering, see below. CG(r) is often referred to 
as translational order correlation or simply 'translational 
order'. Being exact to order e, the average can be raised 
to the exponent as if u(r) was Gau ssian distributedE3. 
We thus can obtain from Eq. (4.8) the result for the 
asymptotic scaling regime, 



CG(r)(x.gGirW(x2 + z2)- 
x(x2 +zf)"2(i+i/'=) 

with exponent 

77G(A^)=A*(«;)(aG)2 
and the geometrical factor 



ig(«) 



(4.11) 



(4.12) 



ga = exp 



l + K V 2 



1 - h 



Zt 



(4.13) 



Our main result on the translational order of the flux line 
lattice consists in the decay of order with a nonuniversal 
exponent 77G(t). Its dependency on k is obtained from 
Eq. (4.12) and is shown in Fig. ^for one of the smallest 
reciprocal lattice vectors G = G p with |Go| — 47r/(aV3)- 
The full decay exponent in Eq. (4.11) is modified by ad- 
ditional K-dependent factors due to the anisotropic elas- 
ticity leading to different decay of contributions from 
transversal and longitudinal modes. The nonuniversal 
behaviour is at variance with the variational ansatz re- 
sults and conjectures in Refs. |lj,|l5|. 

Let us consider the two limiting cases from above: For 
z the result simplifies to 



Cg(x) cx (|x|/ia) 



xexp^A*(K)(aG)2i^((iG)2-i 



, (4.14) 



which again coincides modulo the difference in the co- 
efficent with the result found by GD for this limit. For 
K = 1, the translational order decay is isotropic according 
to 

CG(r')ocLr^'V'r'"'^'^ 

with the rescaled coordinate r' defined above. This 
isotropic limit can be used to demonstrate clearly that 
the triangular flux line lattice considered here and the 
scalar model do not belong to the same universality 
class. Whereas we obtain ridl) — 1.14e for the tri- 
angular lattice, the result in the RG approach for the 
scalar model - which corresponds to a square lattice - is 
,^g(1) = 7rV9e = l.lOe. 




0.2 0.4 0.5 0.8 

FIG. 10. Dependency of the exponent t/Gq on the ratio 
a = caa/cu of elastic constants. 

Another order parameter of interest is the positiorml 
glass correlation function suggested by spin glass theorytZI 



^PG(G,r)-|(*G(r)vI/J,(0))| 



(4.15) 



It measures the thermal fluctuations of the flux lines 
around their disordered ground state. We calculate 
S'pG (G,r) within our framework from the hamiltonian 
renormalized up to scale |r| = K~^e^ . First order pertur- 
bation theory in the pinning energy gives the corrections 
to the mere thermal result to first order in e. We get 



m>l 



Spg{G,v) ^ 5°G(G,r) <! 1 + e ^ c™ ( ^ 



2m 



(4.16) 



where S%q{G, r) denotes the correlation function for the 
pure system with thermal fluctuations only and c„i are 
numerical coefficients. It is finite for r ^ oo and reduced 
with respect to unity merely by the standard Debye- 
Waller factor. As the order e corrections decay they 
can surely not compensate the leading constant term and 
make up for a more than powerlaw decay of the whole 
correlation function. This provides signature of a posi- 
tional glass to order e. 

Widely discussed, however not completely resolved is 
the question, if there exists a phase coherent vortex glass 
state in impure type-II superconductors. A finite value 
of the correlation function 



(4.17) 



for large |r| is ptnppsed to identify such a phase coher- 
ent vortex glasa^j'tJ. Here ^'(r) is the Ginzburg-Landau 
complex order parameter. It can be decomposed in am- 
plitude and phasefactor, 5* = |vl/|e'['^°+*'^l , with ground- 
state phase 00 and phase fluctuations Scf). In the London 
limit the amplitude is constant outside the vortices and 
CvG becomes 



Cvoir) = |*|2|(e«[50(r)-50(O)]^|2 



(4.18) 



Since the phase fluctuations are topologically constrained 
by the vortex positions, the distortions of the vortex lat- 
tice can destroy phase coherence. More concrete, phase 
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fluctuations and vortex displacements are related b; 



D. Experimental implications 



2lT 



(V_L X u(r)) 



(4.19) 



This relation allows for the calculation of Cvg{'^) in the 
framework of the elastic description of the flux line lat- 
tice. We use perturbation theory for Tidis with fluctu- 
ations on scales below |r| having renormalized this per- 
turbation. Then corrections to the mere thermal average 
(. . . )o are given to first order in e by 



CvG (r) = e 
1 



- p-(W(r)-50(O)]^>o 



1 



— ^R*^l [cosh(2i^(r,f))~4cos: 

G 



Wr,f))]}. 

(4.20) 



Here R. 



G 



exp(-GQG/3(wQ(0)u^(0))o) 



e are the 



reduced Four ier co efficients of the random energy corre- 
lator of Eq. ( |2.12| ) and F(r, f ) = r(r - f ) - r(-f ) with 
r(r) = (Gu(r)(5(/)(0))o. The exponential factor of -Rq is 
finite in d > 2. For the disorder induced correction in 
Eq. (4.20) we focus on the isotropic limit with C44 = cgg 



and get 

r(r) = 



T GyX — GxU 
47ra^C44 



l__j^(Ar)), (4.21) 



where Ji(x) is the Besselfunction of ffi'st kind. Upon 
expansion of the cosh-terms, a careful investigation of the 
behaviour for large |r| of the remaining integrals gives to 
order e 



Cvg{v) ^ 
l + e 



,-([50(r)-50(O)]^>o 



diT"^ In(rA) +^2 + 



m>0 



(4.22) 



The di,Ci are again nu meri cal coefficients, yet different 
from the ones in Eq. ( 4.16 ). The exponential factor of 
Eq. (4.20) is determined by mere thermal fluctuations, 
which read 



4^^'^(q)[l-cos(qr)] 
q ^ 



(4.23) 



Corrections to order e can thus not compete against the 
exponential decay of Cyo (r) originating from strong ther- 
mal fluctuations. Thus we conclude that to order e there 
is no phase coherent vortex glass. Whether this result 
is valid to higher orders in e and thus in 3D remains 
however unclear within the present analysis. Dorsey et 
al.cJ indeed found a vortex glass transition in 6 — e di- 
mensions starting from a Ginzburg-Landau Hamiltonian. 
If this transition exists down to three dimensions is not 
clear. 



There exists a considerable number of recent experi- 
mental studies of the structure of flux line lattices in high 
temperature superconductors. Neutron diffraction stud- 
ies provide strong evidence. £fML.the proposed Bragg glass 
with quasi-longrange orderQoLJ. More recently, magnetic 
decoration studies on BSCCO showed hugfi-.dislocation 
free regions containing up to 10^ flux lincaHj. The last 
result strongly supports our elastic description, which ne- 
glects dislocations. 

Comparison with neutron diffraction experiments 
shows the preferred role of the translational order cor- 
relation Cc,{y). The diverging dependency of the cross 
section in a neutron scattering experiment is given by 



CT(k) oc (pk P-k) oc S'(k), 



(4.24) 



where k — (k^, kz) is the difference between in- and out- 
going wavevectors and t he Fourier transform of the 
flux line density in Eq. ( |2.5| ). The structure factor S'(k) 
is the Fourier transform of the density-density correla- 
tion. Close to a reciprocal lattice vector, 

kz=qz, k± = G + q±, q± < G, 

the structure factor becomes 

5(G + q_L,g^)- [ ^e''?^^+*'i^^(e^(G+qx)(u(x,^)-u(o,o))) 

Jz X 

~ J d3re*'^^+"i-^''GG(r). (4.25) 

The scattered intensity is thus described by the Fourier 
transform of th e tra nslational order correlation GG(r) 
defined in Eq. (4.1C). The Fourier transformation can 
be done numerically for general k but one can see easily 
that the power law decay of Gq (r) is slow enough to effect 
powerlaw divergence of S{G + q±, (7z) for small qj^. For 
the limiting cases cn = cqq, k = 1 and cn ^ ceei n — Q, 
of which the latter is realized physically close to , the 
structure factor reads 



^(G + qj 



,qz)(x [cii + —ql 
C66 



(-3+»7g(k))/2 



(4.26) 



As a consequence of the quasi-long-range order in the 
asymptotic regime beyond the scale La, the structure 
factor diverges for small (q^, Qz) leading to Bragg peaks. 
Giamarchi and Le Doussal thereafter named the weakly 
disordered flux line lattice a Bragg-glass. 

In experiments, Bragg peaks and their disappearence 
for higher magnetic fields due to a melting of thejilux line 
lattice have in fact been observed in BSCCOErB. How- 
ever, quantitative details of the divergence cannot be re- 
solved putting our prediction of nonuniversality of "qa 
beyond the precision of today's experimental diffraction 
techniques. But there exist recent experimental results 
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which are based on a more microscopic way to deter- 
mine the structure of flux line lattices. Kim et al. use a 
scanning electron microscope to obtain a spatial map of 
the flux line displacements itself over a region containing 
about 5-10^ linesEa. Due to the small coherence length 
of « 20A in BSCCO we have ^ < a so that only the RM 
and A regime can be observed in principle. Experimental 
results for the mean squared relative displacement of lines 
are available in the RM regime, whereas the asymptotic 
regime is beyond the experimental limit. Kim et al. mea- 
sure in this regime a roughness exponent of C_rjv/ — 0.22. 
This result is in reasonable agreement with our first order 
e expansion result of Crm ~ 0.18 in 3D. 

V. DISCUSSION AND CONCLUSION 

For the first time, the three dimensional Abrikosov lat- 
tice in the presence of point disorder has been treated 
within a renormalization group procedure including all 
of the elastic modes. Triangular symmetry complicates 
the calculation of large scale 3D correlation functions 
in their anisotropic dependency on position arguments. 
These technical difficulties could be overcome and we can 
write correlations in their full anisotropy. More impor- 
tant however, translational quasi-long-range order in im- 
pure type-11 superconductors is described by a - contrary 
to previous claims - nonuniversal power-like decay of the 
order parameter correlations. In particular, the decay- 
exponent rjc depends (if weakly) on the ratio k = Cqq/ch 
of the elastic constants. 

In isotropic superconductors at low temperatures, 
where flux lines interact via central forces, one has < 
K < 1/3. K ~ 1/3 for A < a, i.e. for fields close to 
i/ci, and K — > for iJ ^ Hc2- For most of the field 
region k « (Pq/IGttX^B,. Thus, an increase of the ex- 
ternal field from Hc^ to should result in an increase 
of 7/G and a decrease of (rm- Numerically, the effect is 
small, since ?7go ranges from 1.145 to 1.159 and Crm from 
0.1745 to 0.1763 in this K-range. Thus it will probably be 
hard to detect this effect. However, we are likely to have 
suppressed some of the nonuniversality in our calculation 
when we extended the z-direction to 2 dimensions. In the 
propagator the effect of variations in cn and cge is thus 
reduced by the more heavily weighted C44q2-term. The 
behaviour of 2D disordered triangular lattices supports 
this point of view. Here, Carpentier and Le Doussal ob- 
tain nonuniversal large scale behaviour that is much more 
pronounced than in our higher dimensional casa23. The 
small quantity of the calculated effect should thus not 
mislead to underestimate its qualitative relevance. 

Different from usual RG approaches we observed the 
full flow of renormalization of the the interaction. We 
thus gained in addition to the fixed point dominated large 
scale behaviour information on the intermediate length 
scales. We find a crossover of the structural correlation 
functions from a Larkin-regime, where perturbation the- 



ory applies, to the random manifold regime and eventu- 
ally to the asymptotic Bragg glass regime. Out of one 
Hamiltonian this confirms most clearly the physical pic- 
ture that had originated over many years from different 
approximation strategies to this prominent physical sys- 
tem. For the random manifold regime, where fluxlines 
explore many minima in the energy landscape but do 
not yet compete against each other, we could extract 
the roughness exponent numerically from the flow. This 
is valuable in itself as it can be compared only to an 
estimated exponent from scaling arguments. Moreover 
it also shows dependency on elastic constants, i.e., is 
nonuniversal. 

Dislocations are excluded in our description. Whereas 
the stability of the Bragg-glass phase against these topo- 
logical defects had not always been commonly agreed 
upon, today the position of the Bragg-glass in the phase 
diagram is well established. In the H—T plane at not too 
large fields, weak disorder - the definition of which can 
be cast into a Lindemann criterion - has the Abrikosov 
phase of the pure phase diagram become the Bragg-glass. 
It is bounded by a first order transition melting line with 
negative slope, just like the Abrikosov phase in the pure 
system. This line ends at a critical point. For larger 
fields and small temperatures the Bragg-glass is bounded 
by the prol| ifC | ratioji . of .. v ortices that destroy quasi-long- 
range ordeiillHij'E^liSEilElJ. The exact nature of the large 
field state - it may be called a dislocated vortex glass - 
is not clearly understood, nor is the transition (or mere 
crossing) into it from the Bragg-glass. Berker set up a 
criterion for the shift p£ a first order transition to a second 
order one by disordeio. For a sharp domain boundary at 
the coexistence point he finds a pure first order transition 
to be stable against disorder fluctuations above 2D. This 
is applicable to the melting line where a jump in the FL 
density, which plays the role of a scalar order parameter, 
assures the sharp boundary and a first ordep-transition 
is observed both in pure and impure sampled. 

We would like to thank very much H. E. Brandt, M. 
Kardar and S. Scheidl for valuable discussions. Financial 
support was obtained through Deutsche Forschungsge- 
mcinschaft under grant No. EM70/1-3 (T.E.). 



13 



APPENDIX A: REPLICA PINNING ENERGY 

We want to write the replica pinning energy in com- 
pact form starting from the pinning energy expression 
(2.10). The G ^ 0-terms are the relevant ones, as we 
have argued. 



G#Oa=l 



2T t 



E 



X > e 

GG 



iG[x-u"(r)]+jG'[x'-u''(r')] 



(Al) 



We focus on the x, x'-dependency and suppress the z- 
index of the displacement fields for shorter notation. 



GG' •'^^ 



/X iG[x-u"(x)]+iG'[x'-u''(x')] _ 



X — X je 



/ g-iG[u»(x)-u''(x)]^gi[G+G'][x-u''(x)] 



X / A5(x)e^°'['^-^"''(''+^)-"'(''»l 

J X 

/■ g-iGKlxj-u-'lx)] 

1^ J X 

(a(G)+ E e'[G+G'l[--"'WlA(G'; 



G 

X 



G'5^-G 



— iG[u^(x)-u^(x)] 



(A2) 



A(k) is the Fourier transform of A^(x), it is nonzero over 
a region of size ■ For the first approximation the slow 
variation of [u''(x + x) — u''(x)] with x is used. In the 
region |x| < ^, which is the one contributing to the in- 
tegral, its variation is negligable versus x. In the last 
step, rapidly oscillating terms are neglected versus the 
constant A(G) as they are both multiplied by the slowly 
oscillating e~**^[" i^)] and integrated over. 

This gives the replica pinning energy as stated in equa- 
tion (|2.11|) 



(r)] 



^-^E / ^K(r)-u^(r)). (A3) 



In this version of the replica pinning energy, the effec- 
tive disorder correlator i?(u) is invariant under the full 
symmetry of the triangular lattice. This is because pin- 
ning energy depends only on the positions X + ux(z) 



of the vortices. Redistributing the fluxlines to refer- 
ence positions, i.e., relabelling, may thus not show. Tidis 
in Eq. (2.4) is consistently invariant under ux(2) — 
ux+y(z) + Y with Y a constant lattice vect or. This 
reads u(x) u(x) + Y after substitution (2.8) and ex- 



plains the discrete translational invariance of the Fourier 
sum i?(u). 

The point group symmetries also arise from the invari- 
ance of disorder energy under a change of the 'original' 
positions of the fluxlines in the pure system. Let us again 
resume the distinct notation for the field before (u(x)) 
and after (u(x)) substitution (|2.S|). 



u(x) -> D6oou(x) 



(A4) 



is not an exact symmetry of the pinning energy, it is 
rather the transformation 



I.e., 



u(x) D6oou(x) + (1 - D6oo)x, 



X - u(x) Deoo (x - u(x)) 



(A5) 



that leaves Tidis unchanged, as can be seen in ( |2.7| ) (the 
rotation is absorbed in the reciprocal lattice, that is 
mapped onto itself as Dgo" is a lattice symmetry). For 
the local term G = — G' that we keep for the final com- 
pact form o f th e pinning energy, symmetry ( |A5| ) becomes 
symmetry (A4). It reads for {x, u(x)} 



u(x) ^ u(D6olx) + (Djo - l)x 



(A6) 



The disorder energy in Eq. (2.4) is transformed as 



E fdz^-^Vd^s (X + u(X),z)^ 
E / dz'^-'Vd^s (DeVoX + u (DgolX) , z) 

X 

and is invariant as the lattice can be r elabelled X — > 
DgooX. Our starting Hamiltonian ( ^.11 ) for the RG is 
written in the quantities with the tilde. Correlations cal- 
culated above are thus also expressed in {x, u}. For any 
roughness with exponent smaller than one however, cor- 
relations of u coincide with the ones of u. This justifies 
our dropping of the tilde in most of the treatment above. 



APPENDIX B: INTEGRALS FOR THE 
EQUATION OF FLOW 



In this appendix the coefficients M of equation (3.4) 
are calculated. 

M'-p.-fS ^ ^l-i r g^p{q)g^si~ci) 

J q 



with 



14 



{PapPjS 



al3 
'apPjs) 

The projectors P^'-^ are given after equation (^). Take 
(a/3), (7(5) as column- and row- indices of a 4 x 4 array. 
The indices shall run through {xx,xy,yx,yy). Then 



fp'apP^sig'Y -I 



/ 3 1 

110 

110 

V 1 3 




Here the relations 

»> 



,4 
V_ 

9i 



r2 % Q 



> 



r2 % 



4 

y 



^2 Ixly 

q yj_ 

2 2 



= r gi\{qt+H'ii+^t)= r g 

have been used. Similarily, one obtains 



/; 



plffP^s (g^) 



1 



and 



plpp^sg'-g^ 



/a 1 
Olio 
Olio 

V 1 3 



/ 2 6 

0-2-2 
0-2-2 
V 6 2 





grgL- 



This gives M"'^'^"' as stated in K^) 



( 3/i + 2/2 



M 



with 











/1-2/2 /1-2/2 
/1-2/2 /1-2/2 

V /i + 6/2 




h ^ di-^ /> {g^r + gi) h ^ di-^ /> grgL 

These integrals are evaluated exactly in d = 4 to 
11-1-K, 1 InK 



^ Stt^ C44C66 
with K — cge/cii. 



87r^C44Cii K — 1 



APPENDIX C: FOURIER TRANSFORM OF THE 
ANISOTROPIC PROPAGATOR 

In this appendix the displacement cor relation functio ns 
shall be calculated from the propagator (MQ(q)u;3(q')) in 
Eq. (4.2). The Fourier transformation reads 



Ba/3(r) 



2/ (uc<(q)M/5(q'))(l - cosqr). 



We integrate in d = 4 with the z-direction extended to a 
2D subspace. The sx-displacement correlations then are 



B3;a;(x, Z) = 

= C Jd\ (1 - cos(q^z + q_Lx)) 



i 



1x f2 I % ^2 



= C y dij {1 — cos(q2z) cos(qj_x) + sin(qzz) sin(qxx)} 



'ix ^2 I ^2 
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= C Jd%y Aq^qi_gl y d0 [1 - COs(q,z) 

X cos((7j_a; cos -|- q±_y sin </)) 

+ sin(q^z) sin(gj_a; cos (j) + q±y sin </))] cos^ </> 

-t-f/T-analogue} 

with C = 2A*(27r)~^. The sin(q2z)-term does not con- 
tribute as it is antisymmetric in q^ . So we have 

Sj;a;(x, Z) = 

= Ctt Jd^Qz Jdqj_q^^gl 1 - cos(q^z) 

2 ,2 

x(jo(g±|x|) + ^ J2(g±|x|) 



1 — cos(q2z) 



2 2 

x ( Jo((7±|x|) + ^ ~/ J2(g±|x| 



with J„ the Bessel function of the n-th kind. We used 

„2 



d(f> cos{qj^x COS (f) + q±y sin (j)) COS (f> 







2/2 -x^ 

= Jo((7±|x|) + TT — J2((?_l|x|) 



X 



and its analogue for the t/T-term. With q^ written in 
spherical coordinates as well, the -integration is done 
easily using standard tables, 

Bj:3;(x, z) = 

27r2C Jdq^qj_ J dq^q^ {gl{qj_, q^) [1 - jQ{q^z) 
/ 2 _ 2 

X f Jo(g±|x|)+ ^ ■/2(g±|x|) 

+ gT{q^.qz) [1 - JoiQzz) 

/ 2 2 

X (^Jo(g±|x|)+ ^/ ■/2(g±|x|) 
1 



C44 J Icii 



Z;i4'i(z;qj_) 



X I Jo(g^|x|) + ^ / ■/2(g±|x|) 
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+- 



C66 



- ztKi{ztqi_) 



2 2 

X (Jo(g±|x|) + ~/ J2(g^|x|) 



Kn is the modified Bessel function of the n-th kind. We 
have once more rescaled the coordinates according to zi — 
^■\/cii/c44, Zt — z ^ C66 / C44 . The third terms in both 
the longitudinal and transverse part is not pathological 
{Ki{q) ^ q~^, J2{q) ^ q^ for q ^ 0) and can easily be 
integrated. Then 



-Bx2;(x, z) 



ziKi{ziq±)Jo{q±\K\) 




(CI) 



To separate the cancelling divergences from non- 
diverging 'geometric terms' we introduce the function 



/(7)- / dJ --7'i^i(7''Z)Jo('z) 
Jo 19 



which yields the integral in Eq. ( |Cl[ ) for 7^ = ^/|x|. This 
function can be written as 



/(7) = /(O) 



with 



/(O) 



/•l^l^"' 1 Ixl 
/ dg-(l- Jo(g))=lnU+o(i) 



dqqKoit^q)Mq) 



1 + mj,0^)Ko{t'i^) - t^mMmK.it^m 



2f 



2 |x| - 



1 + 

Integrating the first term of f'{t) gives 

1 



whereas the latter two yield vanishing integrals for large 
Ixl. We now have 



/(7=y^) = lnJ^ + iln l+I^I- 



Oil) 



1.. 



Oil) 



With the definitions k = cqq/cu, h(t) = < ^ In(l-l-t^) and 
considering the rescaling for the equation of flow ( [3.^ ), 
( |3.7|) , as to write in terms of the the dimensionless A* 
from the numerical solution of the RG equations, we fi- 
nally have 



,(x,z)-w,(0,0)]2) = 



-Kin 



In 



La 



La 
2 



La 



2x2 



Zl 
La 

N 

Zl 



All nonconverging terms for z, |x| — s- 00 have been re- 
tained. The j/y-correlations are obtain ed b y k — > 
{x,y) — !■ (y, —x) as explained in Section II B . The mixed 
xy-correlations are calculated similarily and read 

Ba;.y(x,z) = 



X^ C44C66 
.*„2 



,(0)]K(r)-w,(0)]) 



Zt 



Zt 



A*a^ xy 



K x^ 



Zt 
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